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Abstract 

A relationship between the preexponent of the rate constant and the distribution 
over activation barrier energies for enzymatic/protein reactions is revealed. We con- 
sider an enzyme solution is an ensemble of individual molecules with different values 
of the activation barrier energy described by the distribution. From solvent viscos- 
ity effect on the preexponent we derive the integral equation for the distribution 
and find its approximate solution. Our approach enables us to attain a twofold 
goal. Firstly it yields a simple interpretation of solvent viscosity dependence for 
enzymatic/protein reactions that requires neither a modification of the Kramers' 
theory nor that of the Stokes law. Secondly our approach enables us to deduce the 
form of the distribution over activation barrier energies. The obtained function has 
a familiar bell-shaped form and is in qualitative agreement with results of single 
enzyme kinetics measurements. General formalism is exemplified by the analysis of 
literature experimental data. 

Key words: enzyme catalysis, solvent viscosity, Kramers' theory, single enzyme 
kinetics. 



1 Introduction 



The idea of conformational heterogeneity (multiple conformational substates) 
[I], [2], [3], [1] received a new impetus in the last fifteen years with introducing 
the concepts of static [5], quasistatic |6j and dynamic disorder [7], [8], [9], 
PH] . The results of these studies suggest that there is a broad distribution 
over the values of the catalytic reaction rate constant k2 for an ensemble 
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of enzymes and the rate of a single enzyme strongly fluctuates in time [5], 
[Z]; [6], [8], [9], [To]. The reason is that in different conformational substates 
that an enzyme can probe it has different catalytic efficiency. However it is 
widely believed that both static and dynamic disorder of reaction rates are 
essentially indistinguishable in ensemble averaged experiments [7], i.e., the 
distribution can not manifest itself in ensemble averaged kinetic measurements 
and is noticeable only in single molecule kinetics |5], |6], |8], |9], [10]. The main 
reason for the appearance of the distribution is believed to be the fluctuations 
of the value of the electric field strength in the active sites of the enzymes from 
the ensemble [TT]. The constant electric fields in enzyme active sites created 
be preorganized dipoles of protein structure are considered as the key factor 
for the stabilization of the transition state and as a result for the values of 
the activation barrier energies [12], [T3], [H], [T5], [16]. In other words they 
determine enzyme catalytic efficiency [12], [13], [E], [E], [IS]- That is why the 
variability of the electric field strength in enzyme active sites has to lead to a 
broad distribution over the values of activation barrier energy, i.e., to that of 
the catalytic reaction rate constant. 

In our opinion the point of view that the distribution over the values of ac- 
tivation barrier energy can not manifest itself in ensemble averaged kinetic 
measurements is inconsistent. Indeed ensemble studies by Frauenfelder and 
coworkers on rebinding of CO to hemeproteins upon photodissociation long 
ago revealed dispersed kinetics resulted from heterogeneity [1], [17]. The fluc- 
tuations of rate constants (for which recently the term "dynamic disorder" 
was coined [7], [H], [H], [ID]) were also inferred from ensemble studies [I]. One 
can conclude that if the distribution over activation barrier energies exists it 
must affect ensemble kinetics somehow. The distribution acquires a status of 
a unique characteristic of the system and is of primary interest. The aim of 
the present paper is to show that the existence of the distribution is closely 
related with the effective preexponent of the enzymatic reaction rate constant. 
The main premise of the paper is that the existence of the distribution is suf- 
ficient for the appearance of the experimentally observable dependence of the 
enzymatic reaction rate constant on solvent viscosity that manifests itself in 
the preexponent. Thus we persist to attain a twofold goal. On the one hand 
we show that the existence of the distribution enables us to interpret experi- 
mental data for the effect of solvent viscosity on enzymatic reactions. On the 
other hand these data enable us to deduce the form of the distribution that 
can be in principle compared with that obtained from single enzyme kinetic 
measurements. Regretfully this comparison can be only indirect and quali- 
tative because for explored by now enzymes exhibiting the solvent viscosity 
dependence any single enzyme kinetics measurements are impossible. The lat- 
ter are feasible only for enzymes possessing a unique fluorescent active group 
such as, e.g., famous flavin adenine dinucleotide in the cholesterol oxidase \7]. 
However for such enzymes data on solvent viscosity dependence of the reaction 
rate constant are not available. 
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Viscosity dependence of enzymatic and protein (ligand binding/rebinding) 
reactions is known for a long time [H], [19], [20], [21], [22], [23], [21], [25], [26], 
[23, [2B], t29j, [30], [S], [32], [33], [IT], [31], [35], [3S]. For such reactions the 
functional dependence of the reaction rate constant for the rate limiting stage 
k on solvent viscosity 1] has the form 

k oc -a 

iv/Vo) 

where t]o is the viscosity of pure solvent (for water rjo = 1 cP at room tempera- 
ture) and < /3 < 1 (usually (3 ^ 0.4-^0.8). This dependence is experimentally 
verified in the range of variation of solvent viscosity by two orders of magni- 
tude T] < 100 cP. Similar dependence also takes place for folding of proteins 
(see [37], [38], [39] and refs. therein) and at formation of protein structure [ID]. 
However we will not touch upon these processes in the present paper. 

The famous transition state theory in application to rate constants of en- 
zymatic reaction deals predominantly with their free energies (see extensive 
review of Truhlar and coauthors [H] and refs. therein) but fails to shed light 
on the solvent viscosity dependence of their preexponents. The main tool to 
describe the solvent viscosity effect on a reaction rate constant is the high fric- 
tion limit (also called strong damping or overdamped regime) of the Kramers' 
theory [12] or its modifications [13], [11] combined with the Stokes law for 
the friction coefficient. There are many approaches to interpret the effect of 
solvent viscosity on enzymatic rate constant [19], [22], [IS], [30], [31], [33], [16], 
[17] . A brief survey of these approaches was presented in our previous paper 
|47] . The model developed there required neither modification of the Kramers' 
model nor that of the Stokes law. The main premise of our approach was in 
the fact that a realistic enzyme solution is actually an ensemble of individual 
molecules with different characteristics (conditions for the movement of the 
system along the reaction coordinate). It was shown that the experimentally 
observed dependence ([1]) can be obtained if we take into account heterogene- 
ity of conditions in the ensemble. The effective reaction rate constant was 
obtained by averaging of individual Kramers' rate constants over the distribu- 
tion. The aim of [17] was to show that the idea of heterogeneity enables one 
to resolve the problem of solvent viscosity effect on enzymatic reactions in a 
conceptually much more simple way than modification of either the Kramers' 
theory or that of the Stokes law employed previously for interpretation of the 
effect under consideration. 

In that regard we stress that the idea about the importance of nonhomogeneity 
of protein solution for all aspects concerning viscosity dependent effects was 
introduced earlier in the papers [31], [31], [36]. In these papers a useful notion 
of local microscopic viscosity was shown to provide quantitative description 
of experimental data on translational and rotational diffusion of proteins. At 
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this approach the role of hydrodynamic interactions is highhghted. Making 
use of the notion of hydrodynamic radius the authors of [31], [36] succeeded 
in taking into account not only the molecular weight of cosolvent molecules 
but also their shape and size and obtained quantitative description of transla- 
tional and rotational protein dynamics in the presence of polymeric cosolvent. 
However taking into account the nonhomogeneity of the bulk solution is not 
sufficient for interpretation of experimental data on solvent viscosity depen- 
dence for enzymatic and protein reactions because for the latter the protein 
structural fluctuations are of utmost importance [31]. The analysis similar to 
that performed for rotational diffusion is not possible in this case because 
the size of that part of the protein responsible for the fluctuating motion is 
not known [3l]. For enzymatic and protein reactions only experimental data 
for the dependence of fractional exponents on cosolvent molecular weight are 
available [3U] . 

In the present paper we retain the main idea of the model [17] but suggest its 
different technical realization. It seems regretful but inevitable that the way 
by which the distribution can be introduced in the ensemble of enzymes is not 
unique. In [17] we made use of the distribution over the weight with which the 
contribution from solvent viscosity enters into the viscosity for the movement 
of the system along the reaction coordinate. Thus in fact we carried out the 
averaging of the preexponent in the effective reaction rate constant over the 
preexponents of Kramers' rates for individual samples from the ensemble of 
enzymes. In the present paper we explore another option and show that the 
peculiarities of enzyme kinetics enable us to make use of the distribution over 
activation barrier energies. As was mentioned above the existence of the latter 
is convincingly confirmed by experimental data from single enzyme kinetics 
measurements. Thus in the present case we carry out averaging of the preex- 
ponent in the effective reaction rate constant over the exponents of Kramers' 
rates for individual samples from the ensemble of enzymes. The merit of this 
approach is that the concept of the distribution over activation barrier en- 
ergies causes no doubts in the community of the researches. That is why its 
application to the problem of solvent viscosity effect on enzymatic reactions 
seems to provide intuitively more simple and acceptable physical picture. The 
cost for this physical simplicity is somewhat larger mathematical complexity 
of the present model compared with the previous one. Of course the final 
choice between these options will be done on the basis of their theoretical 
non-contradictability and capacity for description of experimental data. 

In in the present model as well as in our previous model [17] the experimental 
dependence of the fractional exponent /3 on the characteristics of cosolvent 
molecules is incorporated as input information. As was mentioned above for 
enzymatic and protein reactions only experimental data for the dependence 
of fractional exponents on cosolvent molecular weight are available [30]. This 
fact imposes a natural limitation on both models. The present model likewise 
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that of |17] is able to take into account only cosolvent molecular weight but 
not the shape and size of cosolvent molecules. 

The paper is organized as follows. In Sec. 2 the discrete averaging of the rate 
constant for the ensemble of enzymes is discussed. In Sec. 3 the continuous ver- 
sion of such averaging within the framework of Kramers' theory is considered. 
In Sec. 4 the obtained integral equation for the distribution over activation 
barrier energies is analyzed. In Sec. 5 the general formalism is applied to the 
analysis of experimental data. In Sec. 6 the results are discussed and the con- 
clusions are summarized. 



2 Discrete averaging of the rate constant for enzymatic reactions 



Our primary aim is to derive the averaging procedure for the reaction rate 
constant. First we consider the value of the reaction rate constant obtainable 
in the experiment. The initial reaction rate for the Michaelis-Menten scheme 

(MM) E + S ^ ^ ES P + E is given by the well known expression [18] 

y ^d[P] _ h [Et] [S] 

dt Km + [S] ^ ' 



where [Et] is the total substrate concentration and the Michaelis constant is 

A-„ = >tl±h (2) 

ki 



The reaction rate constant is measured at the excess of the substrate 

[S] » Km (3) 

Thus one obtains 

V^h [Et] (4) 



Here k2 is the effective reaction rate constant that is the average over the 
ensemble of enzymes. Thus in the experiment one actually measures the value 
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Now let us derive the explicit expression for the value of k^^^ ■ For the sake of 
convenience of notation in this Sec. we assume that the ensemble of enzymes 
consists of a discrete set of pools with different values of the reaction rate 

constants k\. That is we consider the kinetic scheme E ^ S ^ ^ ES^ P + 
E where i = 1,2,3,.... For the quasi-stationary concentration of the i— th 
enzyme-substrate complex {d [ES^] /dt = 0) we have the usual expression 



ES' 



ki [E] [S] 



(6) 



The total enzyme concentration is given by the expression 

[Et] = [E] + Y^ [ESi 



(7) 



From (E]) and ([7]) we obtain the quasi-stationary concentration of the i—th 
enzyme-substrate complex 



ES' 



h [S] [E 



T 



{k., + ki){l + k, [S] 1/ {k^i + ki) } 



Substituting ([H]) into the expression for the initial reaction rate 
d[P] 



V 



dt 



E ^2 [es' 



(9) 



we obtain 



V = k, [S] [Et] E 



+ t^) {l + k, [S] Ej 1/ {k-i + ki) } 



(10) 



At the excess of the substrate we have 



ki [5] E 1/(^-1 + ^0 »1 



:iii 



Under this requirement we obtain 
V ^ [Et] E ^2^. 



(12) 
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where we have denoted 

1 



(A;-i + A;|)E,1/(A;-i + A;| 



(13) 



It is obvious that the values 5i represent normahzed fractions, i.e., 

Y.^^ = l (14) 

i 

Comparing f ll2p with (E]) we finally obtain 

kl" = Y.m (15) 



The latter is the discrete version of the averaging of the reaction rate constant 
for enzymatic reactions. 



Continuous averaging for enzymatic reactions within the frame- 
work of Kramers' theory 



The main tool to describe the viscosity dependence of a reaction rate constant 
is the high friction limit (also called strong damping or overdamped regime) of 
the Kramers' theory |l2]. In it the reaction is conceived as a diffusion process 
of a particle with some effective mass along a reaction coordinate over some 
potential surface. The friction coefficient for the particle is supposed to 
obey the Stokes law /i = Qt^It] where I is the characteristic linear size of the 
particle and rj is th viscosity for the movement of the system along the reaction 
coordinate. The famous Kramers' formula for the high friction limit is 



kx = exp 



(16) 



where Ua and Uf, are characteristic frequencies of the potential surface at the 
bottom and at the top of the well respectively, e is the activation barrier height, 
ks is the Boltzman constant and T is the temperature. 

We denote 

" = 12^ 
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and rewrite the Kramers' formula as follows 

kx = — ; — exp (—ae) (19) 



On the other hand as was discussed in Introduction for enzymatic reactions 
the experiment yields 

k^ff = ^ exp (-aE) (20) 



Here < /3 < 1 is the empirical parameter while b and E are effective values 
for corresponding parameters. As was stated in Introduction the main premise 
of the present paper is to relate effective reaction rate constant k2^^ with the 
existence of some distribution p(e) over the barrier heights e in the ensemble 
of enzymes. In other words we assume that for an individual enzyme (sample) 
from the ensemble the Kramers' formula flTIJl) with a definite value of e takes 
place. Then the effective rate constant (120|) results from the averaging over 
the ensemble with the distribution p{e). Thus the ensemble is characterized 
by the distribution p(e) over the values of the energy e. 

The continuous version of the discrete averaging considered in the previous 
Sec. is obtained as follows: 

kl — > k2(e) = — ^ exp (-ae) (21) 
V/Vo 

E 1/ (^^-1 + ki) ^ Jde f\ ^ (22) 

T J k^i + a/ {r]/vo) exp (-ae) 

oo 

E^i = l^ Jdep{e) = l (23) 



CXJ 

E = J2 ^E = j de p(e)e (24) 
» 

oo 

k'^ff = E ki,6, kl^f = jde p(e)fc2(e) (25) 



Taking into account (!20|) and (fT3|) we obtain from the latter expression the 
equation for the unknown function of the distribution p(e) 



de 



P(e) 



1 + (k-i/a) (ri/rio) exp (ae) 



de 



P(e) 



k_i + a/ {r]/rio) exp (-ae) 



(26) 



We denote 

k-i T] 



a r]o 



(27) 



Then (!26l) can be cast in the form 

1-/3 oo 



tin 



de 







p(e) exp (ae) 
1 + cexp (ae) 



de 



1 + c exp (ae) 



(2^ 



Taking into account the obvious identity 
exp (ae) 1 1 



1 + cexp (ae) 



1 + cexp {ae) ^ 



(29) 



and the righthand side of 0231) we obtain the integral equation for the distri- 
bution p(e) 



de 



P(e) 



1 + (fc-i/a) (r//r/o) exp (ae) 1 + (A;_i/6) (r^/r^o)'^ exp {aE) 



(30) 



The righthand sides of fl23l) and (|2^ yield two equations 

oo 

/ de p(e) = 1 



oo 

E = J de p(e)e 



(31) 
(32) 



for finding two effective parameters b and E that can be directly compared 
with experimentally observable values (note that the parameter a is expressed 
via molecular parameters by f[T7|) and is considered as a given one) . Thus the 
equations (!30!) . (I3T1) and (!32|) represent a closed system of equations for our 
problem. It is worthy to stress that the equation (130|) is strict and exact as 
long as the requirement of the excess of the substrate 



h[S]Jde- 



+ a/ ('7/'7o)exp(-ae) 



» 1 



(33) 
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is satisfied. 



4 Analysis of the integral equation ( 1301) 



We denote 



k_i exp {aE) / a 



(34) 



As we will see later in practice the requirement 
q » 1 



(35) 



is satisfied. In this range the following solution of the integral equation (!30!) 
can be guessed 



asin(7r/3) 
2^ 



1 



+ 



cos(7r/3) + cosh[a/3e — Ing] 



- J2 ( — ] exp{-[l-(l + n)/3]ae} 



(36) 



where [x] in X)'^' means the integer part of x. The latter distribution can be 
verified by direct substitution into (1301) and numerical integration. It yields 
excellent coincidence of the lefthand side of the equation with the function in 
its righthand side for a wide range of physically reasonable parameters defined 
by the requirement (l35l) . In Fig. 1 and Fig. 2 the distribution ( 136|) is depicted 
at different combinations between the values of the parameters (3 and q. 

Within the range of validity of (!36|) outlined by the requirement (!35|) we obtain 
from (|3T|) 



ctg 



71(3 1 



sin(7r/3) 

2g(i-/3), 



o 



1 + q cos(7r/3) 
q sin(7r/3) 



(37) 



From fl32l) we obtain 

TTaE(3'^ = ln(g) arctan 



sin M A r 
^ , +27r/3 n2 -4L — + 

1 + cos 7r/3 y ^ V 2 ; 



L e 
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— arcsm 



1 + cosh (ln(g)) cos(7r/3) 
cos(7r/3) + cosh (ln(g)) 



71(3 sin(7r/3) 

2g(l-/3) 



+ 



(38) 



where 



^ = arctan 



, fH(l)\ 
tanh tan 

V 2 ; 




(39) 



and L{x) is the Lobachevskii function L{x) = — J dt In (cost). Recalhng that 



q = q{b,E) (see (IMIl ) we conclude that equations (137|) and (!38l) are a system 
of two equations for two unknown parameters b and E. 

There are two ways to analyze the equations obtained. The first way can be 
called direct (or purely theoretical) one. At this approach we assume the pa- 
rameter a given by (ITTI) as a known one (i.e., we assume Ua, uJb and / to be 
known values), calculate theoretical parameters h and E from (1371) and (!38|) . 
identify these parameters with effective values of corresponding parameters 
known from the experiment and finally compare the theoretical values ob- 
tained with experimental values. Parameter q = q{b, E) is eliminated as an 
independent value at such approach. In practice this way is hampered by the 
fact that the values Ua, and / are generally unknown and besides it is rather 
difficult to solve the system of equations fl37|) and fl38l) . The second way can be 
called inverse (or experimentally motivated) one. At this approach we identify 
parameters b and E with effective values extracted from the experiment, con- 
sider them as given values while consider the parameters a and q as unknown 
ones. Equation (1371) at large values of q transforms into approximate identity. 
Thus we conclude that within the range of the validity of our approach out- 
lined by the requirement (1351) the normalization of the distribution is always 
approximately satisfied. Then equation (1351) gives the relationship between 
the parameters q and E (direct dependence E = E{q) or implicit dependence 
q = q{E)). In Fig. 3 this relationship is depicted at different values of the pa- 
rameter p. At given values of E and /3 we can define from fl38|) the value of 
the parameter q. From flMl) we obtain that the parameter a is unambiguously 
determined by q at given values of /3, 6 and E 



bq 



exp (aE) 



(40) 



In practice it is more convenient to proceed along the second way. Thus the fol- 
lowing strategy is suggested by the present approach. We consider the empiri- 
cal fractional exponent f3, the effective preexponent parameter b, the effective 
activation energy E and the reaction rate constant k-i as given experimental 
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values. Then we solve the implicit equation ( 138|) to find the value of the param- 
eter g as a function of the values f3 and E. If the obtained value of q satisfies 
the requirement (153]) then our distribution (15^ satisfies the normalization re- 
quirement (1311) . Knowing /3, b, E, q and k_i we obtain from fl4Up the value of 
the unknown parameter a. The latter is directly related with molecular pa- 
rameters of the system via relationship f|T7|) and thus provides information of 
considerable interest. From consideration of f|T9l) and fl20l) we intuitively antic- 
ipate that the value of the parameter a should not differ significantly from the 
value of the parameter b. Then the approximate coincidence of the obtained 
value of a with the value of b would mean that our model works satisfactorily 
well and our approach yields reasonable interpretation experimental data. In 
the next Sec. we exemplify this strategy by the analysis of data for oxygen 
escape from hemerythin presented in [30j . 



5 Analysis of the experimental data from [30] 



There are very few enzymatic / protein reactions for which sufficiently complete 
set of experimental data can be found in the literature. In fact oxygen escape 
from hemerythin explored in [30] seems to be a unique example. Even for this 
protein "glycerol is the only cosolvent for which extensive data are available" 
[30] . For oxygen escape from hemerythin with viscosity varied by glycerol 
kesc = 4 X 109 (r//77o)"°-'^exp [-H,sc/ {RT)] where H^sc = 13 kJ ■ mol-' 
at the temperature of 278 K [30j. Thus we have the following values of the 
parameters: 6 = 4 x 10^, (3 = 0.54 and E = 13 kJ ■ mol~^ at the temperature 
of 278 K The latter means that aE = E/ikeT) ^ 5.63. From Fig. 3 
we obtain for such values of aE and (3 that log q ~ 1.2, i.e., q ~ 15.85. 
This value of the parameter q satisfies the requirement (!35|) though with some 
extension (see Discussion), i.e., our distribution (!36|) satisfies the normalization 
requirement (13T|) . The distribution (l36l) for these values of the parameters is 
depicted in Fig. 4. In Fig. 5 both hands of the integral equation ( 15U]) are plotted 
with the distribution f l36p and the parameters obtained. 

For oxygen escape from hemerythin the role of k_i is played by the rate 
constant kint (s^^) = lO^exp [—Hint/ {RT)] where Hint = 4: kJ ■ mol'^ at the 
temperature of 278 K [30]. Thus we have k^i ^ 1.77- 10^ s^^. Substituting all 
parameters into (HOl) we finally obtain a ~ 2.0 ■ 10^. This value is sufficiently 
close to the value of the parameter 6 = 4 x 10^ that provides strong evidence 
in favor of our approach. We conclude that for oxygen escape from hemerythin 
the combination of molecular parameters is 

^au^^ ~ 2.0-109 



127r2/r/o 
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The latter estimate enables one to make assumptions about the shape of the 
potential energy along the reaction coordinate (characterized by the frequen- 
cies Ua and Ub) or the linear size / of the the effective particle used in the 
Kramers' theory. 



6 Discussion 

We reveal the relationship between the preexponent of the effective reaction 
rate constant and the distribution over activation barrier energies for enzy- 
matic/protein reactions. This relationship arises from the averaging of the 
reaction rate constant in enzyme solution over the ensemble of individual en- 
zyme molecules described in Sec. 3. In the mathematical form it is formalized 
in the integral equation (|30|1 that possess considerable generality and strict- 
ness as long as the requirement of the excess of the substrate fl33|) is satisfied. 
Solution of this integral equation yields the normalized distribution over the 
values of activation barrier energy fl5^ . The typical behavior of the obtained 
distribution over activation barrier energies is depicted in Fig. 1 and Fig. 2. 
They show that depending on the combinations of the parameters (3 and q 
two cases are possible within the range of validity of our approach g >> 1: 
1. the distribution becomes wider with the decrease of the parameter /3 while 
the position of its maximum retains its value unchanged (see Fig. 1) and 2. 
the distribution becomes wider with the decrease of the parameter /3 while the 
position of its maximum is shifted to higher values (see Fig. 2). 

The analysis of experimental data within the framework of the present ap- 
proach requires combined information on both viscosity and temperature de- 
pendence of the reaction rate constant. That is from the experimentally ob- 
servable value of the effective reaction rate constant 



we need to extract three parameters: 6, /3 and E. Besides the value of the 
reaction rate constant A;_i for enzymatic reaction or its analog kmt for protein 
reactions (see previous Sec.) is necessary for the analysis. All these parameters 
are considered as input information for the present model. The theoretically 
calculated value of the parameter a must approximately coincide (or at least be 
commensurable) with the given value b extracted from experimental data. The 
latter requirement is the criterion for the validity of the present approach. The 
example of experimental data for oxygen escape from hemerythin with viscos- 
ity varied by glycerol from the paper [30j analyzed in the previous Sec. testifies 
that this criterion can in fact be satisfied for realistic situations. Regretfully 
these data seem to be unique in regard of their completeness and sufficiency 
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for our approach. New experimental measurements are highly desirable in this 
field of science especially concerning enzymatic rather than protein reactions. 

The present approach retains the main idea of our previous model |17] and 
its capacity for interpretation of experimental data on solvent viscosity de- 
pendence for enzymatic reactions. However technical implementation of the 
present model is different than that of ^7j. Here we make use of the distribu- 
tion over activation barrier energies that has considerable conceptual benefits. 
The latter distribution is not merely an abstract notion. The progress in sin- 
gle enzyme kinetics measurements has done this distribution to be actually 
a directly observable experimental function. The distribution over activation 
barrier energies obtained in the present paper has a familiar bell-shaped form. 
In contrast the distribution over the weight with which the contribution from 
solvent viscosity enters into the viscosity for the movement of the system along 
the reaction coordinate employed in |17] looks somewhat unusual. The latter 
has divergencies (although integrable, i.e., the distribution is still normalized) 
at both ends of the range. Thus the distribution obtained in the present paper 
intuitively seems to be more simple and comprehensible. For this simplification 
in physical picture we have to pay by increased mathematical complexity of the 
present model compared with that from [17|. The integral equation obtained 
in the latter paper had exact solution. For the integral equation derived in the 
present paper we have been able to obtain only an approximate solution. How- 
ever this solution yields excellent accuracy within the physically interesting 
range of the parameters outlined by the requirement fl35p as has been verified 
by numerical integration. We can formulate the range of the validity of our 
distribution f p6|) as follows: the latter is the more satisfactory the better the 
requirement fl35|) g >> 1 is fulfilled. In the case of oxygen escape from hemery- 
thin for which we have rather small value q ~ 15.85 and the requirement fl35l) 
is satisfied with some extension. However from Fig. 5 we see that even in this 
case we obtain satisfactorily good approximation for the solution of the inte- 
gral equation fl5U]) . The small value of q for this protein reaction is a result 
of rather low effective activation barrier energy Hesc = 13 kJ ■ mol~^ at the 
temperature of 278 K [30j that yields E/ {ksT) ^ 5.63. For most enzymatic 
reaction effective activation barrier energies are much higher so that the ratio 
oi E / {ksT) attains typical values 15 ^ 20 at room temperature. From Fig. 
3 we see that for such cases the parameter q attains very large values up to 
10^ -T- 10^. Thus we conclude that for these reactions the requirement fl35|l is 
fulfilled very well. As a result we conclude that for most enzymatic reactions 
our distribution f p6|) provides approximation to the solution of the integral 
equation f l5U]) with excellent accuracy. We anticipate that for such reactions 
our formalism yields very good quantitative description of experimental data. 

The obtained bell-shaped form of the distribution over activation barrier ener- 
gies is in qualitative agreement with that of the distribution over the reaction 
rate constant obtained in single enzyme experiments [TD]. However the search 
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of any quantitative correspondence seems to be meaningless and premature. 
The direct comparison of the distribution obtained from our analysis based 
on solvent viscosity dependence for enzymatic reactions with those extracted 
from single enzyme kinetics measurements is hampered by the fact that these 
two fields have been developed independently and no objects of mutual interest 
have been explored by now. It seems highly desirable to explore an enzyme 
for which on the one hand single enzyme kinetics measurements would be 
feasible and on the other hand sufficiently extensive experimental data for 
solvent viscosity effect on its reaction rate constant would be available. Data 
on such object would provide the possibility for the direct and crucial experi- 
mental verification of the predictions of the present model. In that regard the 
investigations of solvent viscosity dependence for cholesterol oxidase (which 
is the work horse for single enzyme kinetics measurements [7]) and that for 
/3-galactosidase (for which in [TU] the distribution over the reaction rate con- 
stant is obtained) seem to be of primary interest. The experimental data on 
solvent viscosity dependence for these enzymes of the same completeness as 
those for oxygen escape from hemerythin would be invaluable for quantitative 
verification of the present model. 

We conclude that there is the inherent relationship between the distribution 
over activation barrier energies and the preexponent of the effective reaction 
rate constant for the solution of enzymes. Our model yields simple interpre- 
tation and the quantitative description of the available experimental data on 
solvent viscosity dependence for enzymatic/protein reactions. The approach is 
in conceptual coherence with the modern trend stimulated by single enzyme 
kinetics. 
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~ yS=0.5; q=6 




2 4 6 8 10 12 14 



Fig. 1. Distribution p(e) over the activation barrier energies e (eq. (|36p ) at increasing 
values of the fractional power exponent /3 and different values of the parameter q: 
P = 0.5, g = 6 (thick line); /? = 0.6, g = 8; /3 = 0.7, g = 11 ; /? = 0.8, q = 16; 
j3 = 0.9, g = 23 (thin line). 
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— y8=0.5; q=29 




2 4 6 8 10 12 14 



Fig. 2. Distribution p(e) over the activation barrier energies e (eq. (|36p ) at increasing 
values of the fractional power exponent /3 and different values of the parameter q: 
P = 0.5, g = 29 (thick line); /3 = 0.6, g = 41; /3 = 0.7, g = 58 ; /3 = 0.8, q = 82; 
/3 = 0.9, q = 117 (thin line). 
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Fig. 3. Relationship between the parameters q and E (eq. ()38p ) at increasing values 
of the fractional power exponent fi: j3 = 0.4 (thick line); /? = 0.5; /? = 0.6; /? = 0.7; 
j3 = 0.8; (5 = 0.9 (thin line). 
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Fig. 4. Distribution p{e) over the activation barrier energies e (eq. ()36p ) for the reahs- 
tic case of experimental data on oxygen escape from hemerythin from the paper (|30p . 
The values of the parameters are: /3 = 0.54; T = 278 K {a = l/{kBT) « 2.6 • lO^^); 
q = 15.85. 
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left and right hands of eq. (30) 




Fig. 5. Left (thin line) and right (thick Une) hands of eq. (j30p with the distribution 
p{e) over the activation barrier energies e (eq. (f36]) ) for the reahstic case of experi- 
mental data on oxygen escape from hemerythin from the paper (j30p . The values of 
the parameters are: /3 = 0.54; T = 278 K {a = l/iksT) ^ 2.6 • 10^^); q = 15.85. 
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